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Abstract 



Within a quasilocal small sphere, we find the Bel-Robinson tensor B and re- 
cently proposed V have the same gravitational energy-momentum values. One 
may doubt that whether {B, V} are identical? However, we claim {B, V} are 
distinguishable because we find they have different values, in particular, we use 
pseudotensor to evaluate the energy in other small regions such as ellipsoid. 
In addition, as a double check, we find {B, V} have different number of inde- 
pendent components. Finally, we discover the completely traceless property of 
{B, V} implies the conservation of 4-momentum in small region, and vice versa. 



^ ■ 1 Introduction 

In attempts to identify a good physical expression for the local distribution of gravi- 
tational energy-momentum there have been many different approaches which are sim- 
ilar to Einstein's. For example, those of Landau-Lifshitz, Bergmann- Thomson, Papa- 
petrou and Weinberg. Most of them deal with the Einstein field equation G^ = kT^, 
where k := 8nG/c 4 , G^ and T M „ are the Einstein and stress tensors respectively. One 
can define a superpotential with a suitable anti-symmetry U a ^ u = and remove 

a divergence of U a ^ u from G^ v to define the energy-momentum density 

2«* a " := d v Uj* 

O i 

^ \ Using the Einstein equation, we have a total energy-momentum density which satisfies 

d v U a W=2KTa» = 2K{T„>'+t a ' t ), (2) 

o ■ 

C\) ■ and hence, due to antisymmetry of f a , is automatically conserved. Note that T a ^ 

is a tensor, t a ^ is a pseudotensor which does not transform like a tensor under an 
arbitrary coordinate transformation. 

The proposed criteria for testing quasilocal expressions included: (i) limit to good 
weak field values (i.e., linearized gravity), (ii) good asymptotic values both at spatial 
and null infinity. We here emphasize that the criteria for these two are not so restric- 
tive; they actually only test the quasilocal expression to linear order, (iii) positivity 
(i.e., globally) is a strong test but is not easy, (iv) small region: inside of matter the 
quasilocal energy-momentum expression should, by the equivalence principle, reduce 
to the material source terms. Most classical pseudotensors pass this test, (v) small 
region vacuum: positivity for the first non- vanishing parts of the quasilocal expression 
depend on the gravitational field non-linearly and hence this can give a discriminating 
test of the expression, it is quite non-trivial but not impossibly difficult. 

Positive energy should hold not only globally but also quasilocally on any scale. 
However it is generally not at all easy to prove that a particular expression enjoys 
this property. A good test case is the small region limit. This will be our concern 
in this work. Here we consider specifically the pseudotensor expressions. For a small 
region, one can expand the energy-momentum density in Riemann normal coordinate 
(RNC) about the origin: 



Tj(x) = Ta% + d lt r a %af + -dLrJ\ax' t x u + 



Tj\ + d^TJlox^ + ... + tj\ + cy/|o^ + ^tJloxW + ....(3) 



By construction |o and d^ta \o vanish. Consequently, for small x M inside of matter 
the ToP and d^T^ terms dominate (this is a reflection of the equivalence principle). In 
vacuum regions all the T a @ terms vanish, then the lowest order non-vanishing term is 
\d 2 tlv t c P\QX il x v . This is the object of our concern in this work. It turns out that for all 
proposed pseudotensor and quasilocal energy-momentum expressions this fourth rank 
tensor is quadratic in the Riemann (equivalent in vacuum regions to the Weyl) tensor. 
That is why the quadratic curvature expressions become interesting and important 
(i.e., d 2 t a ^ — R a \^a i^V). Normally, the expansion of a pseudotensor expression 
up to second order can only be some linear combination of three tensors {B, S, K } 
or {B, S, W} [1\ which are each certain quadratic expressions in the curvature. 

In the small region vacuum limit, the quasilocal energy- momentum expression, 
which is future pointing and nonspacelike must be proportional to the Bel-Robinson 
'momentum': B /J ,x£ K t tH K (4.2.2 in j2]), where t is timelike unit vector. This is a strong 
test. In the past, B has been considered to be the only tensor which contributes pos- 
itivity in small sphere limit. However, we recently proposed an alternative V (see 
( I24p ) which have the identical Bel- Robinson 'momentum' at the same limit. Unfor- 
tunately, confining a small sphere region, {B, V} cannot be distinguishable which 
implies V is duplicated. Nevertheless, as the basic requirement for the quasilocal 
energy is small 2-surface, there is no restriction it must be a sphere. Examining the 
energy- momentum and angular momentum such as a small box (see Table 1), we find 
they have different values which clarifies {B, V} are distincted. 

2 Quadratic curvature tensors 

Using Taylor series expansion, the metric tensor can be written as 

9*p{x X ) = g*p\ x % + d^g al3 \ x x(x» - 4) + -d^ga^ix* - x%){x v - x v ) + .... (4) 
For simplicity, let Xq = and at the origin in RNC 

9ap\o — Va/3, d^g a/3 \ = 0, (5) 

— 3d^ lu g a /3\o = Raiifiv + Ravfin, ~ ^dyT^ a p |o = R^afiu + R^fiau- (6) 

Three basic tensors [31 H] that commonly occurred in the pseudotensor are: 

B a f3^ u = B( a pn V ) '■= RaXfiaR/3 X u a + RaXucrR^ ^ — „5a(39/iFR 2 , (7) 

o 

SaPiiv = 5'(a/3)(^) = S^vaf} \= Ra^XaRpu^ + RavXaRfifi^ + ~^9a^9 fiuR? > (8) 

A A 3 

Kap^u = i^( a/ 3)( M ;/) = K^ va p := R a \/3 a Rfi if + RaX/3aRu ~ ~9al39^u^ ? (9) 

O 

where R 2 = R pT ^ K R pT ^ K . Using the identity |5] in vacuum R a \ aT Rf3 XaT = \g a pYl 2 , 
some properties of S and K are 

3 

= = Ka^, S a/ 3^ 1 = -(? a ^R = — K a fi^ . (10) 

Unlike B, both S and K are neither totally symmetric nor totally traceless. For 
the non-completely symmetric property (i.e., S a ^ v ^ S afl/3u and K a ^ v ^ K a ^ u in 
general), one can explicitly verify this by using Petrov type D [5]. In particular 

Soon = -7(£n - #11) 7^ 2 (-E ,2 1 - H^) = S0101, (H) 
#oou = l(7E 2 n - 9H 2 n ) + - E 2 U = K 0101 , (12) 



where the evaluation has used the electric part E ab and magnetic part H ab , defined 
in terms of the Weyl tensor [7j as follows: 

E a b '■= C a0 bo, H ab := *C a0b0 , (13) 

where *C a p^ u means its dual. 

It may be worthwhile to emphasize that B has a very good analog with the 
electromagnetic stress tensor T^. In Minkowski coordinates (t,x,y,z): 

T 00 = l - (E a E a + B a B a ) , (14) 
T 0i = 5 il e jab E a B b = {Ex B)\ (15) 
1 r 5 ij {E a E a + B a B a ) - 2 (E i E j + B'B*)] . (16) 



—* — * 

where E and B refer to the electric and magnetic field density. In order to appreciate 
the nice properties of B, we compare some components with S and K 

-Boooo = E ab E ab + H ab H ab , (17) 
Soooo = '2(E ab E ab — H ab H ab ), (18) 
^oooo = -E ab E ab + 3H ab H ab . (19) 

Likewise for the momentum density (i.e., Poynting vector) 

-Boooi = 2£ijkE jd H k d , 5*0001 = 0, K 000i = 2eijkE jd H k d . (20) 

Finally, the stress, 

Booij = 5ij(E ab E ab + H ab H ab ) — 2(E id Ej d + H id Hj d ), (21) 



00ij 



Sij(E ab E ab — H ab H ab ) + 2(E id E d — H id H d ) , (22) 



K 00ij = 5ij(5EabE — ?>H ab H ) — AE id Ej . (23) 

Observing that a linear combination between S and K has the same positive energy 
as B stands: (Aiooo - Soooo ~~ ^oooo) = = (B 00ij - S 00ij - K 00ij )5 v . We define the 
alternative 4th rank tensor [1] as follows 

V :=S + K = B + W t (24) 

where 

3 1 2 

Wapftv := -S a/3 ^ u — -(5g a/3 g Mi , — g a ^g/3u — gaugp^R . (25) 

Both {V, W} have the following properties, let X e {V, W}, then 

X a /3^iu = X( Q( g)(^) = Xftvctfi, Xa/3^ = = X aM /, (26) 

It is known B has the dominant energy property [HI E]: Bap^wfw^w^w 1 ^ > 0, where 
Wi,w 2 ,w 3 , W4 are any future-pointing causal vectors. While V only satisfies the weak 
energy condition and W fulfills none of them. Looking at ( )24"j) . V contains more 
information than B, i.e., V has more non-trivial independent components. For W, it 
does not contribute energy-momentum in small sphere, but does contribute non-zero 
value at the other 2-surface such as ellipsoid or box (see Table 1). We also find that 
it is more convenience using (B + W) to compare B instead of (S + K). 



In our work, we are mainly dealing with 4th rank quadratic curvature tensors. 
There are four independent tensors which form the basis of the space of tensors with 
appropriate symmetries [10J, we may use 



U n 

J -'apfj,v 






p p Act _|_ p p A<7 




•— ^aXpa^^i v T tXaX^a^v fj, 




1 P2 


T 

- 1 - a/Bfiu 





+ -^^R 2 , (27) 

S^^y — -TQapg^v^i (28) 

3 

: #c*/3/^ + o9a09i*>'& 2 > (29) 
o 

(30) 

They are designed to describe the gravitational energy expression based on the pseu- 
dotensor (see fl3Uj) ) and are manifestly symmetric in the last two indices, i.e., M a/3 ^ v = 
Mapfa,,). Then M a p^ u = M^ a ^ v and it also naturally turns out M a p^ v = M^ va ^. Al- 
though there exists some other tensors different from {B, S, K, T} which have these 
symmetries, they are linear combinations of these four, e.g., 

T a ^/3i/ 4" Tocufili = -B a p^ u + —S a i3^ u 4~ ITafi^v. (31) 

In addition, this identity gives another representation for the Bel-Robinson tensor: 
B a f3^u = —-Sap^v + K a /3^ u + -{^gap-g^v — ga^gpu — gaug/3fj,)R 2 - (32) 

2 <5 



3 Energy-momentum tensors of B and V 

3.1 Alternative gravitational energy- momentum tensor V 

Let x^ = (t, x, y, z) and using a RNC Taylor expansion around any point, consider all 
the possible combinations of the small region in vacuum, the total energy-momentum 
density pseudotensor is general expressed as 

Tj = k~ x GJ + {a x B a \ K + a 2 SJz K + a 3 K a ^ K + a 4 f a ^ K )x^x K + 0(Ricci, x) + G(x 3 ), 

(33) 

where a\ to a 4 are constants. Since our concern is the vacuum case, so G a p = = T aj3 . 
Then the first order linear in Ricci terms (9(Ricci,x) vanish. The lowest order non- 
vanishing term is of second order, compared to this in the small region limit we ignore 
the third order terms 0(x 3 ). It should be noted that ToP in ( )33|) is a pseudotensor, but 
in the Taylor expansion on the right hand side the coefficients of the various powers 
of x are tensors. As argued in [10] . d\ v T a (0) must be some linear combination of 4 
tensors, here we use {B, S, K, T}. From now on, the second order term will be kept 
but the others are dropped. There are two physical conditions which can constrain 
the unlimited combinations between {B, S, K, T}: Energy- momentum conservation 
and positive gravitational energy, both considered in the small region vacuum limit. 
First condition: energy-momentum conservation. Consider ([3"3"j) in vacuum 

= dptj = -( ai - 2a 2 + 3a 3 - a^g^R 2 . (34) 
Therefore, the constraint for the conservation of the energy-momentum density is 



0-4 = CLi — 2(3,2 + 3(2-3 . 



(35) 



Although there exists infinite number of combinations which fulfill the above con- 
straint, it has removed one degree of freedom. As each single tensor from {B, S, K, 
T} cannot satisfy the conservation requirement, but a linear combination of them 
can. One can simplify the situation by eliminating T which is absorbed by B, S or 
K. Rewrite (133]) in vacuum 



2 Kjts- 



x^x K 



a>i(B a p£K, + T aj 3£ K ) + a2{S a p£ K — 2T a/3 ^ K ) + a 3 (K a ^ K + 3T a/3 ^ K ) 
= (ciiB a p£ K + a 2 S a p£ K + a 3 K a ^ K )x^x K 

= [(ai + az)B aK n + (a-2 ~ a 3 )S aKK + a 3 W aKK )) x^x K . (36) 

When one considers all the possible expressions for the pseudotensors (some of which 
explicitly included the flat metric), there indeed does appear linear combinations of 
these three tensors [TTJ [T2]. Explicitly, {B, S, K} or {B, S, W}. We prefer the latter 
because W = V — B which does not contribute 4-momentum in small sphere region, 
then B is preserved. 

Second condition: non-negative gravitational energy in small region, i.e., any 
2-surface. For example, we use sphere. For any quantity at t = to as the limiting 
value for radius r where r 2 := x 2 + y 2 + z 2 . The 4-momentum at time t = is 

t p ngivE x = t J x x 3 d x = t ^ij5 3 ^ , (37) 

where / x % x^d 3 x = 4nr 5 /15 when i = j, otherwise zero in spherical coordinate. Thus, 
( 136]) can be rewritten as 

5 

P M = (-E, P) = [{a x + a 3 )B^ + (a 2 - a 3 )S^ 0lJ ] 6 ij , (38) 

where W^oij^ = 0. The energy-momentum values associated with {B, S} are 

B^ = (E 2 ab + H 2 ab , 2e cab E ad H b d ), S^6» = -10(E 2 ab - H 2 ab , 0). (39) 

Here we emphasize that the 4-momentum of B is inside the light cone, — Pq > \P\ > 
0. Observing (|38p . basically we are considering the positive energy within a small 
sphere, B already satisfies this condition and the remaining job is to find {02,03}. 
Equation ( 139]) shows that S^jS^ cannot ensure positivity, since we should allow for 
any magnitude of \E a b\ and \H ab \. The only possibility for (138]) to guarantee positivity 
is when 01+03 > 10|a 2 — a 3 \. However, if we insist the Bel-Robinson 'momentum', 
we only have one choice a 2 = 03. 

3.2 Computing energy- momentum and angular momentum 

Some pseudotensors jlj give the result as a linear combination of {B, V}: Papapetrou, 
2nP a(3 = l(4B a(3 iK - V a ^ K )x^x K . The energy using flE) is 

5 5 

Po = "5iOG (4jBo0iJ " Vmj)5ij = ~lkG Booij5ij > (40) 

where (B 00i j — V 00 ij)8^ = 0. Before we proceed, one might question that perhaps 
{B, V} are equivalent. Does V duplicate? Here we claim they are different since 
Booij 7^ Vooij in general. Elementary, they are constructed from different quadratic 
curvatures: B = B + T and V = S + K + T. In addition, as a double check, we find 
{B, V} have different components which is shown in section 3.3. 



Perfect- 
sphere 


P* = 

jOm _ 


§ *%<^a 5 , r E [0, a], 9 E [0, tt], E [0, 2tt] 
= (0,0,0), ( J 12 , J 13 , J 23 ) = (0,0,0) 


Ellipsoid 


P,= 
JOm _ 


f£ (tVia 2 + tVsfe 2 + t^ssc 2 )^^ x G [-a, a], y E [-b, b], z E [-c, c] 
= (0,0,0), ( J 12 , J 13 , J 23 ) = (0,0,0) 


Hemi- 
sphere 


P,= 
JOm = 

J 12 = 
J23 = 


§ tV^a 5 , r E [0, a], 5 G [0, vr/2], G [0, 2tt] 

- T^(2t°013; 2t°023! t°Qij5^ + t°033) r6 

(^ 023 — t 013) r > <^ = 2l(^ Oijfi 11 + * 033 — 2t 013) r , 
fj^Oy^ + £ 033 — 2t 3 023)^ 6 


Box 


P,= 
JOm _ 


^(^iia 2 + t° M 22& 2 + ^^afcc , x G [-§ , f ], y E [-§, §], 2 G [-§ , f ] 
= (0,0,0), ( J 12 , J 13 , J 23 ) = (0,0,0) 


Cylinder 


P/J- ~ 
JOm = 


ftV/ iJ a 4 /i+i^ 33 ( /l 2_ 3a 2 )a 2^ ; pG [ ,a], y?G [0,2vr], z E [-§, §] 
= (0,0,0), (J 12 ,J 13 ,J 23 ) = (0,0,0) 


Half- 
cylinder 


P,= 
J 01 = 

J 12 = 
j23 = 


fiV/'^+i^^-Sa 2 )^/!, pG [0,a], [0,tt], z G [-§, f] 

4 f /7 5 /i 702 _ liO „3/,3 i 2 (M _i_ 9+0 \ n 5u t03_1 + „3/,3 
Y5 r 012^ ^ — Ys 1 033(3 ft 15v r 011 11 022j^ J ~ g Z 023*2 

Yg^oss^ 3 ^ 3 + 15(^011 + 2t 1 o22 — 2t 2 oi2)a 5 h, J 13 = gt 1 023^ 3 ^' 3 — j^t 3 o\2a 5 h 

Yg(2t 2 023 — ^ 3 033) fl3 ^ 3 — 15(^011 + 2t 3 Q22)fl 5 ^ 



Table 1: Energy- momentum and angular momentum in different small regions 



Usually, the gravitational energy-momentum is integrating in a sphere or ellip- 
soid [13j HU [151 EH] • For angular momentum [3], it is defined as 

J» v := J{xn v \ K -x v e\ K )x^x K d z x. (41) 

We evaluated the energy-momentum and angular momentum in different 2-surfaces 
which is shown in Table 1. The 4-momentum for {B, V} can have different values, 
e.g., box with different lengths. 

Fulfilling the Bel-Robinson 'momentum' requirement using either {B, V} in small 
regions: box allows a = b = c, cylinder and half-cylinder when h 2 = 3a 2 . We observe 
that there is no angular momentum for perfect sphere, ellipsoid, box or cylinder. The 
reason is not the integrand t° ' nijX k 'x 1 'x^ = 0, but the integral / 1° '^jx x l x J d 3 x = 0. 
Conversely, both hemi-sphere and half cylinder have non-vanishing angular momen- 
tum since / t° fl ijX k x l x : 'd 3 x 7= 0. Moreover, we remark for hemi-sphere, if we substitute 
B, J 12 = ^(B 0123 -B 0213 )r 6 = 0. However, if substitute V, J 12 = ^(K,i 23 -K)2i3)r 6 ^ 
in general. Thus, the differences between {B, V} becomes sharp in this case. 

3.3 Counting the independent components of {B, V, W} 

Here we count the components of {B,V,W}. Basically, these three tensors are fourth 
rank, in principle they could have 256 components. However, after considering the 
symmetry properties, they only have a relatively small number of independent com- 
ponents. Counting B has been done, here we claim there is no common term between 
B and W, i.e., {B} f| {W} = {0}. We verify this statement as follows: 

First, we count the independent components of B. In principle, B is completely 
symmetric, by explicit examination it reduces to 35 components. There is a formula 
that directly gives this number. A kih rank totally symmetric tensor in n dimensional 
space has C^ +k ~ l components. For our case C| +4_1 = 35. Since B is completely trace 
free, there are 10 additional constraints which reduce the number of components. 
Therefore, we have left only 25 components for B (another argument see [6]). 

Second, we count V . V does not have the completely symmetric property, but 
fulfill some certain symmetries V a p^ u = V( a p)(p V ) = V^vap as mentioned in ( 126]) . In 
principle, this reduces V to 55 components. However, when we consider the totally 



traceless, the two extra constraints indicated in ( 1261) : V a afJlv = = V a mv . Conse- 
quently we have 55 — 10 — 10 = 35 components for V. 

Finally, we count W. Note that V and W are similar. In principle, refer to 
fl26|) . there should thus be at most 35 components. However, one must also consider 
the extra constraint W a ^^ u ) = 0. This gives 25 additional independent constraints. 
Hence we find 35 — 25 = 10 components for W. 

3.4 Physical meaning of the completely traceless property 

For the quasilocal small region energy-momentum expressions, there are four funda- 
mental quadratic in Weyl curvature tensors. We know that M a/3fJjU = M( a p)(pu) — 
Mfj, va p for all M G {B, S, K, T}. In order to check the completely trace free property 
of the linear combinations of these four, let 



Consider all the combinations, there are two distinct traces because the symmetry: 



8A a ^ au = (ai - 2a 2 + 3a 3 - a 4 )g^ u K 2 , 2A a a ^ u = (a x + a 2 - a^g^R 2 . (43) 
For the completely traceless, the two traces in (113")) have to vanish simultaneously 



Notice that (j4"4"J) is the same as fl35l) which indicates one of the mathematical traceless 
conditions turns out to be one of the physical criteria. Solving ( jUJ) and (|4"5]) . we 
obtain a 2 = a% and this is the 4-momentum requirement denoted in ( I3"8"j) . They are 
1-1 corresponding. From these two results, the completely traceless is not just a 
mathematical property, but covering some physical property. 

4 Conclusion 

Describing the positive gravitational energy, the Bel-Robinson tensor is the best and 
perhaps people may think it is the only one. We recently proposed an alternative V in 
such a way that it shares the same 4-momentum as B does in small sphere limit. But, 
can we separate {B, V} via the different energy value beside this particular sphere 
region? Yes, after examining the value such as ellipsoid, box and cylinder, we conclude 
that {.£>, V} are distinct. In addition, as a double check, we count the independent 
components of {B, V} and indeed they are different. Finally, we discover a necessary 
and sufficient conditions of {B, V}: completely trace free and conservation of future 
pointing non-spacelike Bel- Robinson 'momentum'. 
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